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A coherent coupler is proposed to spin a Bose-Einstein 
condensate composed of ultracold alkali atoms into a vortex 
state (VS). The proposal is based on a Raman transition in- 
duced by two copropagating and a~ polarized Laguerre- 
Gaussian laser beams with diflFerent frequencies. We show 
that the transfer of angular momentum of photons to the 
condensed atoms through a Raman transition leads to a co- 
herent coupling of the ground-state condensate to a rotating 
condensate in a VS. The detection of such a VS is discussed. 



It is well-known that VS's play a central role in charac- 
terizing the superfluid properties of large-size Bose sys- 
tems such as superfluid helium (see, e.g., ef. [Q). Re- 
cently the experimental realizations of Bose-Einstein con- 
densation in trapped ultracold alkali atomic gases |^] 
have generated great interest in studying the superfluid 
aspects of the small-scale trapped Bose gases ||] . To un- 
derstand how the small-scale Bose-Einstein condensate 
is related to a superfluid, it is natural to study the rota- 
tional properties of the condensate and to examine VS's. 
Being different from superfluid helium, the trapped Bose 
gases are not in direct contact to an external container. 
How to rotate such gases and to create VS's is still an 
open question. 

In this letter, we propose to employ a vortex coupler 
to realize such a goal. The principle for the vortex cou- 
pler is illustrated in Fig.l. Two copropagting and a~ 
polarized laser beams along the z direction (gravity's di- 
rection) are used to induce a Raman transition between 
two hyperflne levels of ground-state alkali atoms. Ini- 
tially we assume that a Bose-Einstein condensate is pre- 
pared in the trapping state |— ) = |F = 1, Mp = — 1). We 
consider the time evolution of the condensate after the 
trap is switched off and the Raman laser beams are ap- 
plied. In addition, gravity effects are excluded here since 
we are only interested in the rotational motion of the 
condensates in the x-y plane. To avoid the destructive 
incoherent heating of the condensate due to spontaneous 
decay of excited states the two Raman beams are 
detuned by a frequency A from the optical transition be- 
tween the ground state and the excited state manifolds. 
In this case adiabatic elimination of the excited states 
\j) results in a nonlinear Schrodinger equation describ- 
ing the coherent Raman-type coupling pt] between the 



condensate wave functions 
the magnetic sublevel |— ) and |-|-) 



and corresponding to 
\F = 1,Mf = 1): 
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where ftl^^ denote the positive frequency part of the cor- 
responding Rabi frequencies of the two Raman beams. 
Auj :— UJ2 ~ oJi is the frequency difference between the 
two Raman beams, m the atomic mass, Osc the scat- 
tering length, and TV the total number of atoms in the 
initial condensate, /i denotes the chemical potential 
needed to flx the mean number of atoms. We have nor- 
malized the condensate wave functions by the condition 
/ c?'^a;{|(/)+ P + 10- P} = 1. For the sake of a concise nota- 
tion we have deflned the operator T := ■p^/{2m) for the 
kinetic energy . 

The last terms in Eqs. (|l]) represent the Raman cou- 
pling which coherently transfers the initial condensate in 
the level |— ) to the initially empty level |4-). Evidently, 
if the Raman laser beams are generated by ordinary laser 
sources with Hermite-Gaussian modes, Eqs. (|l|) only de- 
scribe a normal coherent coupler which doesn't cause a 
rotation of the condensate since no angular momentum 
is transferred by the coupling term. Thus, to realize a co- 
herent rotational coupler, the Raman laser beams must 
carry an angular momentum. In view of the recent ex- 
perimental progress in laser physics and atom optics 
a natural choice is to employ Laguerre-Gaussian laser 
beams. For this kind of modes the Rabi frequency Q^j^^ 
for the (7+ polarized beam given by B 
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{I > 1), where r := x + y . A similar expression ap- 
plies to the iT~ polarized beam 5^2^ with ip replaced by 
—if. By analogy between quantum mechanics and parax- 
ial optics, the azimuthal angular dependence of the phase 
factor in Eq. (|^) determines an orbital angular momen- 
tum of Ih per photon in the Raman beams. Here we 



1 



consider the simple case that both Raman beams are 
in the first Laguerre-Gaussian mode {1 = 1) and their 
waist w is much larger than the size of the initial con- 
densate. We then can approximate the above expression 
to fl^i^'^ = \/2i}Q{x + iy)/we'}cp{ikz) and likewise for fl'^^ 
with x + iy replaced by x — iy. Inserting these expressions 
into Eq. (|^) wc obtain 
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We see that, because of the donought-shaped transverse 
profile of the Laguerre-Gaussian mode, the light-induced 
potentials and the Raman coupling terms have approxi- 
mately the spatial structure of a harmonic oscillator po- 
tential V := (l/2)mLi;2^r2, where the effective oscillator 
frequency is defined as ujcS ■= \/nJ(niA)\flo\/w. In order 
to avoid excitations of the condensate after the original 
magnetic trap is switched off, we choose the appropri- 
ate laser intensities so that the effective oscillator fre- 
quency LUeS coincides with the transverse oscillating fre- 
quency of the original magnetic trap. In this sense, the 
donought-shaped Laguerre-Gaussian mode spatially acts 
as a waveguide for condensates where the transvere mo- 
tion of the condensates is confined by the light-induced 
potentials. In addition, one more important feature of 
Eqs. (^, ^ is that the Raman coupling terms now con- 
tain phase factors exp(±2i(/3). 

The physics implicit in Eqs. (|^, ^) is very obvious. As 
sketched in Fig.l, atoms in the condensate corresponding 
to level |— ) absorb photons from the beam and make 
a Raman-transition to the |+) level by emitting photons 
into a~ beam. Each cycle of photon absorption and emis- 
sion in the Raman transition directly leads to the phase 
factors which describe the transfer of an orbital angular 
momentum 2h, originating from two photons in the 
and laser beams, to the condensed atoms. As a re- 
sult, the condensate fraction in the state |+) is set into a 
VS with twice the elementary quantum of circulation by 
the coherent Raman transfer. 

To quantitatively estimate the Raman transfer of the 
non-rotating condensate in level |— ) to the vortex-state 
condensate in level |-t-), we further reduce Eqs. (|[ |4|) by 
assuming separately for each of the two states that the 
spatial part of the condensate varies little during the Ra- 
man coupling. For the present case this approximation 
is justified by two reasons: (1) phonon excitations of the 
condensate 0_ are avoided by choosing the light-induced 
potential V to match the original magnetic trapping po- 
tential, and (2) the spatial dependence of the Raman 
coupling term in Eq. is proportional to which 

exactly leads to a transverse spatial structure matching 



the transverse spatial shape required by a VS with two 
quantum circulations. 

Using this approximation we can make the fol- 
lowing ansatz for the condensate wavefunctions, 
4>-(x,t) = a(t) exp[i(/i/?i — K)t] ipgix) and 4>+{x,t) = 
exp[i(Aw-|-/x/?i— K)t]i/)i,(af), with "^gix) denoting the 
spatial dependence of the nonrotating condensate and "0^ 
representing the VS with two elementary quanta of circu- 
lation. Further, we approximate the spatial parts ipg{x) 
and ^jy (x) of the wavefunctions by that of a trapped ideal 
Bose gas The spatial wavefunctions are then given 
by = exp{-(l/2)[(r/L+)2 + (z/L,)2]}/(^3/4^^^i/2) 
for the ground state and = {x + iy^^g/ {^/2L\) 
for the VS. Here L± and Lz denote the size param- 
eters of the original trap in the x-y plane and the z- 
direction, respectively. The parameter k determines the 
strength of the interatomic interaction and is defined as 
K := TTnascN/{m{2TTf/^L\Lz). Projecting Eqs. §, |) 
on these two states we arrive at a system of two ordinary 
differential equations, 
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where jap gives the occupying fraction of the conden- 
sate in level |— ) and \I3\^ that of the VS in level |+). 
Because of probability conservation, we have the con- 
straint |a|2-|-|/?p = 1. The terms proportional to Woff/v^ 
describe the laser-induced Raman coupling between the 
non-rotating condensate and the VS. The nonlinear con- 
tributions proportional to k have their origin in the inter- 
atomic interaction. As the strength of the nonlinear in- 
teraction depends on the density of the condensate, which 
is different for the non- rotating condensate and the VS, 
the nonlinear terms are of different magnitude for the 
two states. Finally, the expression linear in (3 of Eq. (||) 
includes the energy difference 2fiuJcs between the VS and 
the ground state as well as the energy transfer Til^uj from 
the laser beams to the condensate. 

To solve Eqs. (||, |^) it is convenient to consider the 
population difference f{t) jap — which also is 
the physical quantity we are interested in. By using the 
initial condition /(O) = l,/(0) = (all atoms are ini- 
tially in the |— ) state) one can derive from the system of 
equations ^ the first order equation 

= (1 - /) [^^l^il + /) - ^(1 - /)(19 - e + 9/)^ 
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with £ := 4(Aa; + 2uJcs)Ik- Eq. indicates that f{t) 
can be expressed in term of elliptic functions 0. The 
exact analytical solution is given by 
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dn(z|TO) is one of the Jacobian elliptic functions. The 
parameter a is determined by one of the real roots of 
a 6th order polynomial P{d = Yl^i=o9i^^ ^ where gi are 
2nd order polynomials in the parameters e, k. The other 
parameters h, m, and v can be expressed as a rational 
polynomial of 5th order in the parameter a. The explicit 
expressions are rather lengthy and will therefore not be 
given here. Instead, we will rely on a numerical calcula- 
tion to determine a, b, to, and ly. 

As an example we use the experimental parameters of 
the ®^Rb experiment |^. The number N of condensed 
atoms is 4500, and the laser intensity to be chosen so 
that the effective oscillator frequency Woff is equal to the 
transversal trap frequency = 132 Hz. The trap size 
parameters are given by L± = 2.35/^m and = lAfim. 
The scattering length is taken to be a = 5nm so that the 
nonlinear coupling parameter becomes k = 422 s~^. In 
Fig. 2 we show three solutions corresponding to different 
values of the frequency difference Aw. The dotted line 
corresponds to Aw = 2700 Hz for which the parameters 
in Eq. (|) take the values a = 0.86, b = 0.39, m = 1.07, 
and v = 514 Hz. As will be explained below the dot- 
dashed line (Aw = 1962 Hz, a = 0.72, b = 0.70, to = 
656, I' — 16.4 Hz) is a very peculiar case for which the 
maximal population transfer is reached (for fixed Aw). 
Within our model about one half of the *^Rb atoms could 
be transferred to the |-f ) state. The dashed curve shows 
that for a slightly smaller frequency difference (Aw = 
1850 Hz, a = 0.70, 6 = 0.74, to = -2.97, u = 236 Hz) the 
solution does only perform relatively small oscillations. 

In order to better understand this behaviour we con- 
sider the equation which determines the extremal points 
in the time evolution of the population difference, f — 0- 
According to Eq. (|^) these extremal points are deter- 
mined by the (real) roots of a 4th order polynomial in 
/. This allows us to determine the maximum population 
transfer for a given Aw. In Fig. 3 we show the real roots 
of the polynomial given by the r.h.s. of Eq. (0) as a 
function of Aw. / = 1 is always a root, and since we 
start just with this value (all atoms in the |— ) state) the 
solution can only oscillate between / = 1 and the next 
root for a given Aw. It can be read off from Fig. 3 that 
the maximum population transfer is about 50% and ap- 
pears around Aw = 1962 Hz (solid line in Fig. 2). The 
shape of the oscillations for frequency differences around 
this value is heavily affected by the fact that the poly- 
nomial becomes very small in the middle of the cycle, 
causing a large delay in the time evolution. For slightly 
smaller values of Aw a new root appears relatively close 
to 1. This is responsible for the fact that the maximum 
population transfer is dramatically reduced if Aw passes 
the critical value around 1950 Hz. 

Though a population transfer of 50% is satisfactory it 
would be desirable to transfer almost the whole conden- 
sate into the VS. Fig. 3 indicates that this is only possible 
if we can reach the peak at / = — 1 (corresponding to a 



complete population transfer) . As is shown above this is 
impossible for a fixed Aw if the system starts from / = 1. 
But this restriction may be circumvented by introducing 
a linearly time-dependent frequency difference between 
the two lasers. We have numerically solved Eqs. (||, 
by assuming that the frequency difference varies with 
Aw(t) = 2900 Hz - (3336 Hz )'^t. Obviously this leads 
to almost complete transfer of the condensate to the VS. 
This effect can be explained by noting that the nonlinear 
interaction between the atoms causes a time-dependent 
energy shift between the two states. This energy shift 
is due to the fact that the nonlinear interatomic inter- 
action is proportional to the populations which change 
with time. Thus, after some time the transition will be 
out of resonance even if it was initially at resonance. A 
phase modulation of Aw can compensate for this time- 
dependent energy shift. 

So far we have shown how to create a VS by coherent 
Raman coupling. How to detect the VS remains to be 
discussed. A direct detection of the VS is to observe the 
off-resonance absorption image of the rotational cloud. 
For a vortex state one should expect a bright " hole" in the 
image which accounts for the vortex " core" in the density 
distribution of the VS. Another possibility to detect a 
VS is to observe the Doppler frequency shift due to the 
quantized circular motion of the atoms. Here we focus on 
discussing the absorption image for the circular motion. 

To realize a non-destructive detection scheme we pro- 
pose to employ two counterpropagating weak and a~ 
probe beams in the x-direction and to observe the ab- 
sorption image of the a~ probe beam. The principle for 
the detection is simply based on Raman absorption of the 
beam (if the VS is prepared in the level |+)). The 
Raman absorption coefficient of the a~ probe beam is 
proportional to the probability of the Raman transition 
from 1+) to |— ) which has the form |^ 

r°° lo |2 
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with W y^S^ + |f7+|2|rj_|2/(4A2). n± are the 
Rabi frequencies of the two probe-beams. S — Aw — 
2{fik\/m) — Awd is the effective frequency difference in 
the two probe-beams with Awi, = AiTikL jrn^ex • being 
the Doppler frequency shift induced by circular motion 
of the atoms in the VS. Eq. (||) clearly shows that the 
Raman absorption image directly displays the rotational 
properties of the VS. Fig. 4 shows the image of a VS 
discussed in this paper. We see that the absorption is 
asymmetric around the vortex axis. Such an asymmetry 
directly gives signature for the Doppler effect induced by 
the quantized circular motion of the atoms in the vor- 
tex state. Since the probe beam propagates along 
the x-direction, the strong absorption (black spot in Fig. 
4) corresponds to atoms with velocity antiparallcl to the 
wavevector khfix of the beam. The weak absorption 
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on the opposite side of the vortex Hne appears because 
there the rotating atoms are moving paraUel to the a~ 
beam. 
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Figure captions: 



Fig. 1: Two Laguerre-Gaussian laser beams induce a 
Raman transition and generate a VS in the condensate. 
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Fig. 2: The population difference between the two 
ground states as a function of time for different constant 
values of the frequency difference between the lasers, and 
for a linearly varying frequency difference (solid line). 

Fig. 3: The extremal points of the population difference 
as a function of Aw (in Hz). The time evolution starts 
at / = 1 and moves downwards until it hits another ex- 
tremal point. 

Fig. 4: A typical Raman absorption image of a VS. The 
Doppler shift of the rotating atoms leads to a higher ab- 
sorption if the atom move towards the a~ laser, y and z 
are given in units of the trap size parameter. 
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